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1.0 SUMMAaJH 


^ mont4i Combustor Lloer Durability Analysis Program (NAS3- 
2183o) perKonnGd4)y Pratt- St—Mhitnsy Aircraft ar© — presGnted, T-he objGct-i-v© of 
the program is to- evaluate the use of advanced three-dimensional transient 
^ u noe-linear stress-strain analyses- for modeling the-thermo- 

mechanical response of a. simulated combustor liner specimen. CycUc life pred.- 
iction technology for creep/fatigue interaction is evaluated for a variety of 
state-of-the-art tools for crack initiation and propagation.. The sensitivity 
of the initiation models to systematic changes..- in- the operating conditions_is 
also assessed. 

The- selected component for the program is- a simulated combustion chamber outer 
liner.. Testing under accelerated cyclic thermal loading provided- the- struc- 
tural response and failure- data .for comparison with the analytical simulations. 

Transient and . steady state- three dimensional heat transfer analyses were con- 
ducted using a finite difference solution..aood correlation wiith measured tem- 
perature response was obtained.- The greatest area of uncertainty in the analy- 
sis was m the calculation of the heat transfer coefficients> Modification of 
the mitral values was required to match the steady state temperature distri- 
bution. 


Results of the transient- heat transfer analysis were used to create an incre- 
mental temperature file for the three dimensional non«linear finite element 
analysis of a-segment of the combustor liner. 


The three dimensional non-linear finite element analysis was conducted with 
the MARC computer code. The analysi s- used existing time independent classical 
plasticity theory with a Von Mises yield surface and the combined (iso- 
tropic-kinematic) hardening rule. A constant rate creep model was used to 
account for instantaneous time dependent plasticity effects. Both the plasti- 
city and creep models were calibrated to isothermal Hastelloy X material 
response data. Results of the analysis indicates that this material model does 
not accurately predict the stress-strain response under thermomechanical load- 
ing conditions. The prediction shows continued cyclic hardening while experi- 
mental data suggests that the response is stable after several-loading_cycles. . 


Three existing state-of-the-art crack initiation models were considered for 
predvction of the combustor liner fatigue life. The Strain Range Partitioning 
method oy-erpredicted the observed cracking by more than a factor of 2. This 
cd ( cul dt i on i s sensitive to tho dofinition of the generic fatigue curves aiid 
the prediction of the inelastic strain components. A Pratt & Whitney Aircraft 
(Conmercial Products Di\dsion) developed Certbustor Life Prediction Wethod also 
overpredicted the crack initiation life. This method is similar to Strain 
Range Partitioning but differs in the definition of the inelastic strain com- 
ponents and the basic fatigue curves. In actual design practice, this pro- 
cedure is used along with extensive experimental and field service data to 
provice a realistic estimate of overall service life. The third method, the 
Continuous Damage Approach, underpredicted the observed cracking by more than 

mechanVeaW a'tigue ot Mgjfy rStl'depende'nt al°s 


A linear elastic fracture mechanics 
growth observed in the combustor li 
crack growth rate than observed. 


approach is used to predict the crack 
ner. All calculations predict a slower 


2JL-IJ41ROOUCTION 


The overall -aperafeing cost of the modern- gas turJalne engijie is stgn1f1cantT.¥_ 
affected by the- durability and. efficiency of the-major hot section coif>ponehts» — 
These are the combustor and turbine structures- in the- engi-oe^ during each, 
flight cycle, these components undergo large- tiiermally induced stress and- _ 
strain cycles which include significant -amounts of creep and relaxation-,. 
Primary r-esponsdbilities of the combustor, in the engine cycle, are gas teih»- 
perature level and pattern control-, required for effic.i.ent turbine operation, 
and- exhaust- emission control at the various fligh.t. operating conditions. These 
goals are. accomplished b.y. the precise metering, of- air^ throughout the combustor 
structure. The. high pressune and- high combustion, gas- temperature character- 
istic of this environment require that the combustor liner be coolad-for dura*-- 
b-ility. These- requirements for control-of exit gas temperature, emissions,, and 
metal temperature generate an intense competition for uti-li-zotion of combustor 
airflow.- The more agqresSive performance, efficiency, and emission- goals- set 
for current and. future engines- emphasize the need for development of- durable 
combustor structures which can operate with reduced, levels of cooling air. 
This requires detailed knowledge of the operating environment and the ability 
to accurately predict structural response for these 1 oadi ngs .... 

Current structural technology activities for hot section component design and 
life prediction seek to exploit major advances that have been made in non- 
linear structural analysis capabilities. Large, general purpose finite element 
programs such as the MARC code have been developed expressly for nonlinear 
structural analysis. These roqrams involve sophisticated computational algor- 
ithms and advanced finite element formulations, yet rely on material models 
whose applicability to the hot section component environment is questionable. 

Of primary concern is the response- of materials to cyclic loading involving 
simultaneous creep and plastic behavior. A major need is the development of 
appropriate hot section component structural response data, sufficient to 
evaluate the advanced structural analysis capabilities with emphasis on the 
effectiveness of the material models. 

Adequate prediction of the nonlinear stress (strain) response of the hot sec- 
tion component also requires good prediction of the thermal load history. Such 
predictions must be capable of modeling both spatial and temporal response, 
suitable for life prediction purposes. Creep-fatigue life prediction for hot 
section components is currently calibrated to simplified models for predicting 
the local cyclic hysteresis response to thermal-mechanical loading. Further, 
these models generally lack calibration to well-controlled hot section compon- 
ent fatigue test data. Thus, a second major need is for the evaluation of life 
prediction models for creep-fati que response of hot section components using 
the results of nonlinear stress (strain) analysis and well-controlled compon- 
ent response data. 

The use of advanced nonlinear models for three-dimensional stress (strain) 
analysis of hot section components poses major burdens in terms of engineering 
and computing times. Design of these components requires cost effective, 
shortcut tools which can evaluate the effects of design changes on component 
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3.0—TASK I - HOT SECTIOa-OlMPiXtENT SELECTION 
3,t~.Combustoc L1r>er Rig 

At Pratt & WiJitney Aircraft,, sirtiul at-eci combustor lioers-have been evaluated in 
a specially constructed experimental rig ► Cyclic engi-ne level temperatures and 
thermal gradients are- generate on siroolated’ combustor liners by a 250 kW,. 450 
kHz induction beater used in. conjunction with contr4>lled cooling air temper a- 

ture*and frlow rate, A schemati c~drawi ng of- the r1g_shQwtng the BLlncipal.com- 

ponents-Js shown in Figjicfi_l. 

Cooling air iS supplied to the lower plenum, of the rig from a non-^itiated 
upstream, air. heater at temperatures up to SSS^C (lOOOOF).. The air in- the 
lower, plenum is then, directed through flow- straightener plates to the upper 
flow- divider plate which permits ratio! ng shroud-Side airflow to- hot-side- 
tlow.- Shroud- side air provides the total cooling air suppli'^d to the specimen 
which- is discharged into the hot-side flow annulus and then^ .together w-iih the 
metered hot-Side flow, ejected to ambient.. — 

The specimen is inducti vely- heated, by power suppl ied .f rom a 250 kW, 450kHz 
induction heater. High frequency (450 kHz) heating has been chosen to minimize 
the depth of penetration of the induced current (heat) in the tes^ specimen to 
best simulate the thermal surface loadings (radi ati ve/ conductive) experienced 
by combustor liners in engine operation. . 

The induction coil is supported on aufromework of glass-bonded mica which is 
transparent to the induction field. A 48.3 cm (19 in) diameter quartz cylinder 
(also transparent to the field) is positioned between the coil and the speci- 
men to form the inner boundary of the hot- side flow annulus. The quartz cylin- 
der is. retained in position by a 45.4 kg (100 1b). quertz cover. The test 
specimen selected for this proqrarn simulates a combustor outer liner and 
incorporates five (5) sheet metal Hastelloy-X louvers-. Prior to testing, the 
test section, transition pieces and top shroud cover are assembled as a unit. 
This provides a permanent instrumentation installation and facilitates fre- 
quent inspection during testing.. The assembled louver test specimen and cover 
are shown in Figure 2A. Fine wire chromel-aliiftel thermocouples are- used to 
determine temperature distributions during testing. Typical thermocouple loca- 
tions are shown in Figure 28. 

3.2- Cyclic Thermal Loading and failure Data 

For this program, results of a single cyclic test were used to provide the 
temperature and failure data for the thermal/structural analysis and life P>*6“ 
dictions. All information from the test was taken from the third (center) 
louver of the test specimen. 

The test program consisted of a ninety (90) second thermal cycle in which the 
rig coolii.g air was maintained at a constant temperature and flow rate 
(504OC-2.5 kg/ sec (940OF-5.5 lb. /sec.)) while cycling power from the 
induction heater. A representative transient and steady state temperature dis- 
tribution measured on the center (third) louver is shown in Figure 3. The 
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^ Sweaty (20) second traasi-ent from an isothermal (504PC 
temperature distribution, a for-ty (40) secorKi steadt 
state-portion, and a-cool-^owo., back to the..or4 ginal . i^t^^^ 

^™e?Lunf.SS‘S\*^^ war» a,erage,O0 pmvide a cir- 

A total, of. 1730 cycles were accumulated, on the- test s-oecimen Insoectinn fnr- 

consent and f-luor ope net rant (zyglo) techniques ^started at 
cycW 500 and continued in intervals of 125 er-2^ cycles The arL 
of the smallest detectable crack, that the eddy current method can reliablv 
determine on^tje test s^Je€imen was 0J)254 cm (o:oi0 Tn^ Onfe 

r ck was determ ined ,_crack lengths_were measured using: the zyglo techniq ue. 

At cycle 1000 detectable cracKs were disco-vered at the edae~ of the Ioqvp^ 

orinin<f ihspectton revealed additional crack- initiation 

ongyns at various circumferential locations at tl^e- edqe of the^^^^^^ 

c?Ick®locUior?b^^^^^^ propagation. g??wth mode. /f. typical 

c^ack Wh data for fnuf in.Fiqur.e 4. Representative 

crack locations are shown in Figure 5. No other 

iSon if tho°f observed during the test. Metallurgical exaS- 

frarture surfaces revealed predominantly a tr.ansgra^lar crack 

t^TJrack behi?^^ surfaces prevented further SZiZ lf 


Based on these results, axial lip cracking was 
mechanism and location which would be addressed 
life prediction. 


identified as the key failure 
by the structural analysis and 
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Measured Thermal Response 
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TASK U - HEAT TRANSFER AND STRESS-STRAIN ANALYSIS 


4.1 Heat T»^attsf^r Analysis 

4.1.1 Thermal Analysis Method 

transient three-dirtiensi orval heat tOT 

test specimen was. conducted- to nft»*fiQns of the program. The anal- 

the structural. TwhiSy Aircraft developed thermal analysis 

ysis was- performed using a Pratt &vmiy geometric breakup prepro- 

systar.. The "laj or elements of- this lo^^er into a 

cessor which. prxivtdes_ a . whose properties ar.e associated viiith. a 

series of individual thermal JyJ!® specifies- environ- 

centroidal nodal ^ heat transfer coefficients, along the sur- 

meot at. temperature, . a^lvses^ both the internal proper- 
faces of the model. For toecific hea^ and heat generation.) 

ties at the nodes (thermal. coefficients, environmental temper- 

and the. external boundary and heat- generation) were specified 

atore, fluid temperature,, fluid ° then obtained through a 

as fu^tions of time. The temperati^e at the nod^ equation at each internal 
finite difference solution the^nternal temperature sol- 

point... Surface temperatures with the struc-- 

utions and the local was^maintained, to a 

ic p!"eafcud 

4.1.^ Model Description 
The discretized section of 

(Jimenstora of the ■"oda' ««« that the measured 



Snet r"e7^ .d a^u^n. coo.h. ho,es ,s 

higher than the remainder of the model. 

In the weld Ion, perfect 

gTpTs^Vsled^o^e^^^^^ upper and l.^er sheet nodes over the initial 

Ind final 0.22 cm (0.085 in) in the weld region. 

4.1.3 Thermal Boundary Conditions 
Convection 

There is convective heat loss ^Th^ on^^except^^^^^ at the end 

;lt^"^^J;^rrd^:"r;\hrla^^ spedm/n is connected to up and 
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downs tiv am louvors. In ordor to provide a continuous conduction path witliout 
IwvJmi to model all the- louvers in tlw" Uist specinieH»_J:he corresponding up- 
stream ami downstream nodes are thermally connected. 

Convective heat transfer coefficients for the hot and cold sides, lip cavity 
and cool ini] holes wore calculated w4tb-the channel flow ecpiation: 

N(\nj = .023 Npiv ^'^NpR 

where Nussult number 

NpR = Praodtl fuimber 
''*Ro “ Reynolds number 

The hot side coefficient was bas<sd upon- the coolinp flow- and slot dimensions 
determined by -the distance between the hot side of the specimen and the quartz 
liner. Tiie heat transfer on the cold side of the lip is enhanced by the 
impintieifient flow. An empirically- determined cross- flow impinqement corre-lation 
was used to calculate the heat transfer coefficients tn this region* - 

Radi ation 

In addition to the convective heat transfer, the upper and lower surface nodes 
radiate to the rig walls at the cooling temperature. The shape factor was cal- 
culated from t1ie equation for concentric cylinders. The emissivity of the test 
specimen and outer rig wall was assumed to be 0.8 while tlie quartz inner wall 
was taken as 0.9. No radiation is assumed froili one part of the model to the 
other, such as within the lip cavity. 

Heat generation was assumed to be circumferentially uniform. Therefore, there 
was no heat transfer into or out of the model in a circumferential direction. 
The otily-circumferenti al variation introduced is by the cooling hole geometry. 

Power Input 

Heterminati on of the heat flux on the louver was accomplished by a separate 
experiment. With the cooling air turned off,, the specimen was subjected to a 
rapid step Increase in power from the induction heater, A high speed recorder 
was used to monitor the initial thermocouple respuanse. By assuming no heat 
loss from the specimen at the thermocouple location, either through convection 
or conduction, the rate of temperature rise was related to the heat flux by: 

« - .. d r 

9 - P ^'p^ THl 

whore c., = metal heat capacity (joule/kg 
q = heat flux (WATT/m^) 

T = metal temperature 
V s metal thickness (cm) 

0 T time (sec) 
p I metal density (kg/m-^) 


n 


Thf» measured- heat flux values from the experiment were then us£d t& cal cii-l^e 
fiijy di^rlbutiort alonq thfi Vouvsr dutilug the actual cyclic test-. dis- 
JrtbSiMTs 5to>S In PUiTe r bi. tte' analysis, thn- heat flu*. »» 
introduced to the specimen thcough the hot side surface nodes. This- ® 
approximation of high frequency induction heatijig where the pflW£ai_inBiL 
rapidly with distance into^the^ specimen, 

ThP variation of input power (heat flux) with time for the cyclic fatigue test 

was contmr.ed by a DAT/VTRAK system. This ^ S 

Figure 8.. During the-. transient heat transfer analysis the flux was applied in 

a sfmilar manner. 

4.1.4 Results of Heat Transfer Analysis— 

The heat transfer analysis- of the louver section was completed i^^ 

Tnitiailv a steadv state calcul atTon- was performed using the heat transter 
coefficients and heat flux, val-ues deterifiined by the methods discussed in sec* 
Sm 4 1 3 ThSe results indicated that the predicted temper-atures were cooV 
"IhVn'flie t5,lmocoupU data. The groatest ^Wer^ce between t e 

orediction and the data, at any one location was 360C (65^). this 

analysis the greatest degree of uncertainty was associated with the 
tioi of the heat transfer coefficients. Since 9r®ater accu^cy was requi^^^ 
fnv> 1-hp <tfriictural analysis and life prediction, modification of the initial 
heat transfer coefficient values was made to improve the agreement, ^ 
avfrai iTc^ase of 5^^^ the originally computed values was required to 
MUrth^mMS^ ed temperatures. The ?'®‘Tibution of heat transf^ 

:Sr?reSrc?ll: ‘a^d r^’roulle'STa 

area of the cooling holes.- 

For the transient analysis, the modified l«3V‘"Th/s®^s’]TeaSi‘f «funS° 
ranain constant throughout the transient eye e. This is a reason.ble 
tion since the cooling air temperature and flow rate chained constant our g 
rhfcvclfnq The comparison of the predicted transient temperature response 
a“ tte measure? diT ls shown In Figure .11.. The analysis sccuraUly predicts 
the transient response at all locations within the louver. The 
of the thermal c^le include the maximun temperature gradient occurring early 
in thi S (10 seL) and the inverse temperature gradient between the edge 
nf the louver lip and the weld regions during the cool down portion of the 
cycle. After the initial heatup on cycle 1, all subsequent thermal eye es were 

predicted to be identical. 

The final results of the transient heat transfer solution were used to create 
an incremental temperature file for the structural analysis. 
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4.2 Sti^es^Stral 0 -Analysis 

4*2U_ Non»U near .Str4ict.ura1. Aj9aXy.s4 s -Met hod 

Ttie MARC general purpose finite element program (Reference- 1) was used, to pre- 
dict the response of the louver component, to the thermal- and mechanical loads. 
In evaluatijig the response the program assumes that the- loading history is 
divided into a number, of increnentally applied ioading steps. Each, load- step 
is then sequentially analyzed -as a linear problem using an appropriate stif- 

fness matrix and load vector. 

The plasticity algori thm..i:n MARC uses a tangent stiffness approach, in.v»hich. 
the stiffness matrix, is updated and reassembled for -every plastic loading, 
increment.. An Iterative -scheme is then, used to calculate the- r.esponse of the 
structure- to each individual load Increment. When the- strains used to estimate 
the stiffness matrix in the assembly phase- are- equal,, within- a. user specified- 
tolerance, to the- strains obtained- in the- solution phase, the iteration iS 
terminated and the program, proceeds to the next loading increment. If.con.vergr 
ence is not obtained within a number of cycles specified by the user, the pro- 
gram will terminate with, an appropriate error message. In order to ensure that. 
the solution does not drift away from a state of equilibrium as loading pro- 
gresses, the applied incremental load. vector, AP, is modified by the addition 
of a residual load correction vector, P-I. This residual load vector is ob- 

tained-by taking the sum of the internal reactions, I, and subtracting this 
vector from the total external load vector, P. The residual load vector from 
the (n-l)th increment is then used to augment the applied incremental load 
vector in the nth increment. 

For general use in structures undergoing small deformation, the program^ uses 
the Prandtl-Reuss stress- strain equations of classical plasticity with a 
Von-Mises yield surface. The yield surface may be assumed to expand isotropic- 
ally using the ISOTROPIC option. This option allows the material to harden, 
but due to the equality of the yield condition for tension and- compression, it 
does not allow for a Bauschinger effect. The KINEMATIf option allows the yield 
surface to translate without expansion but predicts a closed hysteresis l^p 
after the first cycle under strain control. The COMBINED hardening option 

allows the yield surface to expand and translate and was considered as the 
best model to represent the cyclic material response. A piecewise linear 

representation of the monotonic stress-strain curves was provided through the 
MARC- user subroutine WKSlP.. 

The nonlinear algorithm employed in the program assumes that the total strain 
may be additively split into elastic, plastic, and creep^ components.. Tne 
plastic strain epl i& determined during the plastic loading phase accord- 
ing to rate-independent plasticity theory. It is asswed that the prior creep 
strain epp does not affect the subsequent plastic response, so that no 

interaction effects are included in the constitutive theory. 
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The creep oetiorv allows time dependent efCects to be ijicTuded by assuming that 
the creep response equations- are constitutive equations whicn are valid—not 
only f’or constant stress histories but for general stress histories. The 
creep strain Is assumed- to. be Independent of~ prior p-lastlc deformation. The 
program algorithm assumes that the nonlinearities due. to the creep constTt^ 
tive relation” may te incorporated into a suitab^le- Toad vector and 
an. Initial strain techn-ique. Ibus» during snvall deformation creep the stif- 
fness matrix remains constant and onJy the Toad, vector- is changed- on an incre- 
mental has 1s^ Since - the stiffness- matrix is- constant, the. reassembly of this 
nratrix Is avol-ded and the program solves for the Incremental displacements by 
matrix back substitution. It has been demonstrated- that this rtesul ts in. t^- 
ical incremental conputer timOs which are about 50% of the— computer time 
required for the assembly and solution required In a plastic loading locre- 
raent. An automatic time stepping scheme Is- Incorporated Into the- solution 
algorithm, to take- advantage- of the dif-fusive- nature of the- creep process. TJie 
description of the material creep response for this analysis was incorporated 
within the MARC user Sibroutlne-CRPLAW. 

Under the high temperatures encountered during the thermal loading cycle the 
material will exhibit time dependent plastic behav.1or - the material ci^eps 
during the plastic loading phase.. In order to properly account for such behav.- 
lor a time dependent plasticity theory is needed.. In default 6f such a theory 
the behavior was approximated in the MARC— program by applying increments of 
creep and plasticity alternately. 

4.2.2 Louver Finite Element Model 

The three dimensional (3D) finite element model of the louver used in the 
structural analysis Is shown In Figure 12. The model contains an included 
angle of 0.577^ which represents the distance between the radial syi^etry 
olartes of two adjaicent cooling holes. This small model size was selected 
because. of the circumferentially uniform temperatures in the louver lip »*egion 
and to minimize computer run time. Aspects of the modeling philosophy,, which 
had-tabe established prior to the complete analysis, are discussed -separately. 

Element Type 

As shown, the model consists of 30 elements. Twenty-one (21) of the elements 
are the 20 node isoparametric element with reduced integration using 8 Sauss 
points (MARC element #57). Since it is known that some reduced integrat on 
point elements are unreliable if distorted shapes are used, the fu 1 27 Inte- 
gration point elatient (MARC element #21) was used around the cooling holes., 
Selection of the reduced integration point eleiflerits was based on initial esti- 
mates of the savings in run time over the full integration point elementr, 
to minimize the required spatial temperature interpolation between 
transfer and structural analysis meshes. 


the 


and 

heat 
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Boufutar^ Condittons. . 

FORMATjnw pernvltted*. Th.ts was. acco<l\pl Ished by use af the TRANS- 

fSre sJs^^.^'Along°''w!e^^^ AB ANlf'^S' effect 

^AB .* ^ CO X 



where; R^3 ori girtal.. radial .coordinate oa.AB 

Rcd orfgimi wdlaJ coordthate. on_C0 

Double Precision 

investioltt^'^tuh'^ problems associated with the small included angle were 

£‘s-er7Ss"i? 

ii:^°^si,ij^s'ir:ra'’drbT‘’ 

the roun“d‘.’^?,-e^rt’o ptSucVa^Vcll^p-taSre^'irot^e^m'tfsoV^^ 

•Ct'^“'®l’'**®5ration point element in double precision was 

a°rh?o7gh1h1irn^s iiH 

tic^^solution indicate a^lmost exa?t“;j^e“^e°^t' ith^^“’fno™ insJi’cli; ^il“l 
Mesh Density 

?he*l“v2 “!' T^'las'VnUiallloStl^^ an optimma mesh along 


15 


cons 1 dered^ -The res-ults of the analyses (shown In E-lgure 13-) Indicate Idttle 
sensittvUy of the predicted hoop ( circumferential^ mechanloal stra.l4» d1.s4n?1- 
butlon to the number' of elements in. the model... Based on tt>ese- results the- 
e1ght..elemeotjaesb-.(mesh #-1) was selected for the 3D -louver- model—. 

Weld Simulation 

In the actual louver, the seam weld extends only over a potd:1on of- the over- 
lapping sheets- (see Figure 6)». During the heat-up portion of the test- cycle 
the hatter leuvec lip wants to expand radially a larger amount than, the colder 
knuckle and cooling hole region. This results analytically in overlapping, of 
the structure In tivis region (an impossibility).. Analyses conducted, using the 
TYING option in.MARC to prevent the overlapping, of these nodes indicated a-13%. 
Increase In. hoop mechanical strain relative. .to the case where- overlapping was_ 
allowed to -occur. 

The phil-osophy used for the 3D model- was to tie- the nodes to prevent overlap- 
ping during the heat-up and steady state portions of the cycle and then use 
the TYING CHANGE option, to untie the nodes during the cool-down portion of the 
cycl e . 

Initial Verification of 3D Model 

Prior to running the complete thermal cycle, the 3D louver model was run with 
Isothermal and pressure loading conditions to verify the construction and 
modeling assumptions. The free thermal (isothermal) test case, with all 
required boundary conditions, transformations, and ties required for the small 
included angle model, predicted essentially zero stress throughout the struc- 
ture.. The second test case was run with a uniform pressure load (2.67 x 10^ 
N/m2 (4 psi)) simulating the actual pressure drop that occurred in the com- 
bustor liner test. Results of this analysis indicated a maximum hoop stress of 
-6.89 MPA (-1 KSI) and are in good agreement with the predicted stresses 
obtained from an axisymnetric louver analysis with the same pressure loading. 
Based on these two cases, the constructi.on of the model was verified.- 

In the final form, the 3D model required 640K double precision words of stor- 
age on the IBM 370/3033 system^ 

Thermal Increment File 

The loading sequence for the- louver analysis consisted of ini tally applying 
the mechanical pressure load followed by the temperature increments as deter- 
mined from- the heat transfer analysis. The thermal loading cycle was repeated 
as required. As discussed in Section 4i,l.l, the heat transfer finite differ^ 
ence map was determined based on locating a temperature solution point as 
close as possible to the finite element integration point. Thus the 4:1 ele- 
ment relationship between the two meshes shown in Figures 6 and 12. The 
results of the transient heat transfer analysis discussed in Section 4.1.1 
were used to define the temperature history throughout the louver model. With- 
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in the MARC program*, the-all^jw^le tamper ature increment s-ijte is estimated 
from the relation: 


AT 


I-' 


la 


where: - A T 

Oy 

E 

a 

F 


= temperature ijicrement (OC) 

= yield stress (MPA) 

» Youngs s M^oduT us (MPA) 

= coefficient of thermal expansion (0C*1) 

= factor based_on in-house experience (0.30 in this 
anal-ysi s) 


This, gui-deline is used to prevent excessi ve strain (stress) changes occurring 
during, a. loading increment possibly effecting the stability of the solution.. 
(Generation of 'the temperature- increment file and required spatial, interpola- 
tion between the heat transfer and finite element meshes was accomplished by 
the general interpolator described in Reference 2. Initial experience for the 
3D model with this interpolation was not successful, particularly in the area 
of the overlapping sheets. The final thermal file required a considerable 
amount of manual modification of the temperature increments to obtain tempera- 
tures at all integration points to within 30C (5°F) of the- desired value 
from the thermal analysis. In the complete analysis, the successive tempera- 
ture increments were obtained from the thermal file with the MARC subroutine 
CREDE. 


4.2.3 Material Model 


Cyclic Plasticity Model 

The combustor liner Sf^cimen, discussed in Section 3.2" was constructed of 
Hastelloy X sheet material. This is a fine grained nickel base alloy which has 
relatively low strength but high ductility at elevated temperatures. It is 
representative of the combustor liner materials used in the gas turbine 
industry. Previous experimental testing conducted,..both internally by P&WA and 
under Reference 3, have shown that tbis material is highly strain rate (time) 
sensitive at elevated temperature, t.e.., creep and plasticity occur simultan- 
eously. The temperatures attained in the cyclic louver test are well into the 
range where time dependent plasticity will comprise a significant portion of 
the structural response. In addition* the differential heating and cooling of 
the louver produces a variable strain rate throughout the loading cycle. The 
general material model available in the MARC code (and most other non-linear 
codes) is based on classical time independent plasticity theory with a Von 
Mises yield criterion. Cyclic plasticity is considered by either the iso- 
tropic, kinematic, or combined hardening rules. For the louver analyses, a 
cyclic plasticity model consistent with an estimated average strain rate v#as 
used together with a creep simulation to account for time dependence at the 
lower strain rates. 
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The development of the HasteUoy X material plasticity model considered a 
tri-linear representation of the- monotonic tensile curves and the COMBINED 
cyclic hardent4>g rule. The criteria for development of the model was based on 
prediction of the cyclic inelas-tic strain range observed during stabilized 
isothermal strain controlled testing at representative strain rates-. In con- 
structing a monotonie curve for each temperature-within the analysis the elas- 
tic modulus, work hardening slope and stress level (at large strains-) was 
chosen- based on previous testing conducted at 0..008 min.“l. The yield point 
was. determined using the tafigency point of the actual stress-strain eurve and- 
by . equating the areas under the experimental and ana.lytlcal curves.- Figure -14 
shows a represent atijffi construction. Using this approach., the monotonic 
stress-strain curves- for temperatures between- 427QC ( 8 OOOF) and 9820C 
(I 8 OOOF) were constructed and incorporated into the MARC WKSLP user sub- 
routine.. The complete- stress- strain-description used in the louver analysis is 
presented in Figure 15, 

A series of isothermal, uniaxial simulations were run to verify the accuracy 
of the representation under monotonic and. cyclic loading. Figures 16A through 
16C show the predicted results at 76QOC (HOQOF), 8710C (I 6 OOOF) and 
9820C (180QOF) for cyclic strain controlled loading between +0.KK strain. 
Shown are the stable 5th cycle MARC. results vs. available cyclic test data for 
two representative strain rates (e = 2.4 X 10"2 MIIi”l and e* 2.4 X 10"3 
MId*l). As shown, the stable stress values and the cyclic plastic strain 
range (width of the loop) are consistent rela.tlve to the test data.- — 

Creep Model 

A creep solution was included in the analysis to model the time dependent 
material response at the slower strain rates. This was accomplished using the 
MARC subroutine CRPLAW and required that the incremental plasticity solution 
be periodically stopped and "creep allowed to occur." 

The final form of the creep simulation incorporated into the CRPLAW subroutine 
is: 



^CR 

= Auf^t 


^CR 

= creep strain 

where: 

0 

= stress (KSL) 


t 

* time (hrs) 


A, n 

= temperature dependent constants 


The simulation of the instantaneous time dependent nature of tne high tempera- 
ture material response required that the constants be determined from short 
time (<1 minute), high stress level (aapplied > 0,5 ayield) HasteUoy X 
creep data. 
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The accuracy of tWs constant rate model was demonstrated by sim«3atlon of 
raonotontc stress re^ajcatlon tests ^ Stress reTaxatien». not creep», is appro- 
priate in judging the-unoddl- s-i4ice». with the thermal loading,, ttie time-depend- 
eot louvef response is predominately stress relaxation. Comparison of these 
results are shown in Figure 17, Good agreement is obtained with- the- model , 
paf-ticulac-ly for the hlghfitJfcamperatures (^8710C (I 6 OOOF)).- _ 

A piecewise interpolation-scheme- was developed for the temperature dependence- 
of the constants. A- and n^ In tte louver analysis, creep was ''allowed'* to occur 
between 1300.QF and 180Q9F. The values of the constants were determined as4 

704?C-871°C 
(130GOF-1600OF) 

871 oe -9270c. 

(1600OF-1700OF) 

9270C-9820C 
(170QOF.-1800OF) 


where T = metal 
4.2.4 Results 

The fatigue critical location of the combustor liner specimen is the edge of 
the louver lip. The observed failure mode is axial crack propagation from the 
edge toward the weld. Results of the structural analysis indicated that this 
location experiences essentially a uniaxial (circumferential/hoop) stress 
field. In the following discussion, the predicted hoop stress/strain compon- 
ents at the element integration point closest to the edge have been used, to 
define the response. 

The predicted results for the first two thermal loading cycles are shown in 
Figure 18.. A representative temperature history used in this analysis for 
three locations in the louver model (end of the lip, weld, and knuckle) is 
shown in Figure 19.. Referring to the letter designations and the temperatures, 
the response at the-end of the lip throughout the two loading cycles can be 
followed, 

A - (0 sec - 5040 c (940OF)) Start of cycle 1 

B - (5 sec - 7320C (135QOF)) Initial yield point 


A = 3.06 X 10-12 /_L_\ n = a.39 X 10-3(T) 

1 1000 1 

A - 9.48.x 10-30 /jJ03.64 g 3 10“g (T) ^31. 49 

\iooof 

A = 2.49 X 10-23 /jl_\ ^5.78 ^ ^ -1.2Z.x 10-2( +25.39- 

llOOO I 


temperature (OF) 
of the 3D Non-Linear Analysis 
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C - (12»5‘ sec. - 9210C (1&9Q0F)) 


D - (6Q..sec...... 954:9^^ 


E - (66 sec --760OC (14000F)) 
F - (78 sec - 5130C (9550F)) 
a' - (90 sec 504OC (940OF)) 


B@ginn4-Rg; of tirst cneep increment. The 
increase- In mechanical strain from B to C 
is produced. ty the increas-ing^ temper^iture 
difference thmugiiouf. the. leaver . The 
reduction-. In- stress level is associated. 
w-ith... the temperature- dependeoGe- of "the 
stress- strain curves used in the. analysis. 

End. of the heating portion of the thermal 
cycle.. Temperature between C and 0 remain. 
essentially constant.. The response 
between. C. and D.. consists of creep/reTaxa-* 
tien and elesti-c unloading. Rationale for 
three CS) creep increments used is dis- 
cussed below. 

Reverse -yield point during cool down. 
Minimun strain point. 

End of cycle 1, beginning of cycle 2. 


Due to the tensile residual stress at the end of- cycle 1» subsequent loading 
in cycle 2 requires more thermal loading increments to produce re- yielding of 
the louver lip (point B').. The rapid temperature rise in the heatup portion of 
the cycle results in this yield point occurring, at 893°C (1640°F) vs 
7320C (1350°F) as in cycle 1. The lower stress level at B'' 159 MPA-(23 
KSI)) reflects the temperature dependence of the stresS-strain curves in the 
analysis. The remaining points in the second cycle (C - F‘) occur at similar 
times as in the first cycle. 


As discussed previously, the analysis was conducted with the time- independent 
plasticity formulation in the MARC code. The stress-strain curves that are 
used, were generated at a strain rate representative of the expected average 
rate experienced during the loading cycle. To account for the actual time 
dependent material response^ a creep solution was included in the thermal 
loading cycle. In attempting to determine the amount of creeps time and at what 
point in the loading. cycle it should be included, the variation in straiorate 
throughout the cycle was considered. Figure 20 shows the predicted variation- 
in strain rate at the end of the louver lip for a portion of the first thermal 
loading cycle. The strain rate was evaluated based on the predicted mechanical 
strain increments and the time associated with the temperature change for that 
incr^ent. The results showL a rapid (>.01 min-1) initial strain-rate, a 
peak value occurring at 6 sec and a dramatic drop in strain rate after 10 
sec. Since the predicted response during most of the first 10 sec of the cycle 
is elastic loading (A - B, A' - B') and not generally rate sensitive, only the 
remaining 50 sec of the heat-up was considered for the creep simulation. The 
50 sec of creep time was arbitrarily divided into 3 segments and applied as 
follows: 
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5 sec (10-15 sec.) 
15 sec. (15- - 30 sec) 
30 sec (30 - 60-sec) 


app^lied at the— 12.5 sec poijit- 
applied at the-20 sec-point 
applied at the 45 sec point — 


As shown-in Figure 18, the three (3) creep so-lutieps. produce botti- a reduction 
in stress- (relaxation.) and. an increase in strain. Elastic unloading 
each, creep sol-U-tioo. The sensitivity of these- results to 
was. investigated by including all of the creep time (5° 

pptnt ia a separate l-oading; cycle-. The results indicated *• change (6%.) 

in the total inelastic (plastic and creep). strain pred4ct6d-.f or the cycle 


The second cycle (and all remaining cycles) required 78 MARG increments as 
follows:— 


60 thermal loading (35 heat-up, 85 cool down) 

14 creep 

3 nuff (no^lSd -to.. equilibri un prior, to creep, and. after the tie change) 

Each loading cycle required 45 minutes of CPU time on an IBW 370/3033 computer. 

The louver model was finally run through a total of six (6) complete- thermal 
loading cycles. The predicted hoop stress and strain, history at the end of the 
lip for the 6 cycles is shown in Figure 21. As- shown,. cycles 2-6 have similar 
stress-strain loops that , are- moving (ratchetting) in the 
direction and show a progressively higher peak tensile stress during 
down portion of each cycle. Details of the sixth cycle response shown in 
Figure 22 and contain the same elastic, plastic and creep regions discussed 
for cycle 2. The total mechanical strain range is defined as the distance 
between points C* and F. The analysis required 4 loading cycles for the total 
strain range to stabilize. This result is shown in Figure 23 where the i|^t]a] 
stSn ranS of 0.34061 (cycle 1) increases to 0.34821 by cycle 4 and^ remains 
constant for cycles- 5 and 6. This smaJl change in the total mechanical strain 
range with the first 4 loading cycles is apparently produced- by a decrease in 
the stiffness of the louver. With the combined hardening cyclic plasticity 
model used in the analysis, a number of loading cycles are required to promice 
a stabilized stress-strain response at each location in the louver. The reduc- 
tion in local stiffness, associated with the change in slope of the ^tj'ess- 
straln curve from the initial to stabilized curve, in combination with the 
pressure and thermal loads on the louver-, produce tha 2 % change in mechaoieal 
strain range. 

The results of the first six (6) cycles also show that the cyclic inelastic 
strain range is decreasing with each cycle. This can be seen in Figure 24 
where the individual inelastic straih components are plotted ys. cycle number. 
The compression going inelastic strain consists of the time- independent plas- 
tic strain generated between points B and C and the creep strain produced 
between points C and 0.. These individual components occur during the heat-up 
portion of each cycle and together comprise the total comprepion going in- 
elastic strain. The tension going inelastic (time- independent) strain occurs 
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between points E amt F during the cool-down portion of each cycle*. As shown, 
the total inelastic strains are not exactly fully reversed which results in 
the rachetting of the stress- strain response*. The change in inelastic- strain, 
range Is consistent with the fact- that the stress range is increasing in each- 
cycle. Frortv Figure 21, the peak tensile stress produced during the cool por- 
tion of each cycle- is increasing. With the higher peak tensile stress, the 
louver lip location now. has a longer elastic range to traverse before yielding 
can occur during heat-uja of. the next' cycle. T-hls requires a larger temperature 
difference (more loading increments) which, reduces the- amount of temperature 
'•load!^ available to produce plastic strain. The high temperature stress levels- 
do not change significantly after cycle 1 because the- cyclic stress-strain 
curves have- har-dened to their final valme,. The predicted creep ctrain per 

cycle remains essentially constant, consistent with the stable stress levels. 

Stress- strain results at other locations throughout the 3D louver model are 
presented in the following Figures. Results from the 6th loading- cycle indi- 
cate: 

1. The maxiituni cyclic strain, range- occurs at the end of the louver lip. 

Figure 25 indicates a rapid drop off in cyclic strain range along the lip. 
These results are consister»t with the observed fatigue failures originating at 
the end of the louver lip in the rig test. 

2. Cyclic inelastic strain in the lip occurs in only the last 3 elements. 

Figure 26 shows the 3D louver model with J<ey elements numbered. Figure 27 is a 
composite stress-strain result for the last 3 lip elements (28, 29, 30). A 

small amount of inelastic strain is observed in element 28 produced by the 
creep solution, otherwise the response is elastic. 

3. Elements closer to the weld region show reduced values of residual tensile 
stress. la Figure 28, element 24 shows a completely compressive stress 
response. 

4. Results in the knuckle region (Figure 29) show: a- peak tensile stress dur- 
ing the heatup portion of the loading cycle and a residual compressive stress 
at the end of the cycle, as expected.- The cyclic strata range at this point is 
approximately 65X of- the cyclic range at the end of the lip. 

4.2.5 Therittomechanical Specimen Test 

As shown in Section 4.2.4, the predicted stress-strain response at the end of 
the louver lip had. not stabilized after 6 thermal loading cycles. The mean 
strain and stress values were increasing and the amount of cyclic inelastic 
strain per cycle was decreasing. Accurate fatigue life prediction for this 
thermomechani cal response requires a detailed knowledge of the history of 
these and other response parameters. Obviously the selection Of the material 
model (Section 4.2.3) determines, to a large extent, the predicted response. 
While the modeling assumptions were calibrated with isothermal test data, the 
actual component response under varying strain rate and temperature conditions 
may be significantly different. Testing of a smooth, uniaxial, strain control- 
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led specimen was conducted to provide material response data- un-der tfe. thermo- 
mechanical conditions for further- model calibration, and to define a stable 
response for tlie Ufe-predIctiooHJortlon of the program. 

P-r-att & Wh.1tney A4r-craft conducts thefTfiomechanical- fatigue testing on- a 
routine basis^ The. experimental system developed to conduct thes« tests Is 
based uporr stand-ard- closed loop Scrvohydraullc test machines using low fre- 
quency (lOkHz) Induction heating and. compressed air for temperature control ort 
the specimen (Ref.. 4). Specimen temperature measurement is provided by radia- 
tion pyrometer and used In conjunction, with Independently computer controlled 
preprogrammed mechani caJ strain- and. temperature histories. The system is cap- 
able oi cycling at a positive or negative mean, .strain, with either stress- or 
strain control and hold times within each cycle. The tubul ar specimen used for 
this type of testing is shown in Figure 30. Axi^l strain is measured, from the 
machined Internal ridges utilizing a linear variable differential transformer 
(LVJ>T) and a quartz internal extensometer. Load and total strain for ttiis sys- 
tem can be controlled to within 155 and temperature within 2°C (40F) 

The predicted circumferential mechanical strain and temperature histories for 
the 6th thermal loading cycle at the edge of the louver lip were used to 
define the specimen test conditions. This- is shown in Figure 31 where the let- 
ter designations are consistent with the points- In the loading cycle described 
In section 4.2.4. Because the region near the edge of the louver lip exper- 
iences a uniaxial stress field, tho stress-strain output from the specimen 
test is considered representative of the actual response producing the fatigue 
fail ure. 

During the first 14 cycles of the test the strain range was inadvertantly set 
larger ( Ae~ .39% vs. = .35%) than the desired value. This resulted in 
the slightly "fatter" stress-strain loops shown in the data. After cycle 14, 
the strain range was corrected. A total of 138 cycles was accumulated on the 
specimen with the "non-linear" strain temperature history. General observa- 
tions-from the test results are: 

1. The stress-strain response stabilized within the first complete cycle 
(cycle #3, cycles 1 and 2 were required for start-up) 

2. No appreciable relaxation of the response was observed -during the 138 
cycl es 

3. The hold time had essenti ally no effect on the shape of the response 

4. Reverse plasticity is observed during the cool down portion of the cycle 
Results from the testing are presented in the following Figures: 

Figure 32 - Comparison of cycle 3 vs. cycle 11 indicates a stable response 

Figure 33 - Cycle 11 vs. cycle 20 showing the effect of larger strain 

range during the first 14 cycles 

Figure 34 - Cycles 60-70-80 indicates repeatibil ity of the response 
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The comiiapison of the experimental results with- the predicted 6th cyc-U 
response shown In Figure In general, the predicted- res p^se— for this 
cyc-Te shews qualitative agreement with the experiment. Further investigation 
of the ability of the-material model to- reproduce the experiJliental- results, was 
conducted with... a one dimensi ona 1 , strain. controlled, analytical simulation ot 
the experiments. This analysis used the. same material model as the 3D louver 
analysis and- was run through the same mechanical stralor temperature history as 
the specimen test. The simplicity of the analysis allowed a larger number" of 
cycles to be econcrolcally run to evaluate material cyclic response. These 
results, In conjunction with-results from previous isothermal analyses, result 
In the foil owning conclusions relative to the material model. 

1. Under Isothermal,, strain controlled cycling, the combined (isotropic- 
kinematic) hardening model will, predict a stable stress- strain- response-. 
This Is demonstrated In Figures 16A, 16B,.16C. 

2 The inelastic strain produced by a creep and/or stress relaxation does 

not Increase the size of the yield surface. This may result in an 

apparent rachetting. of the stress- strain -response. 

3. For the thermcmechanical cycle considered, the predicted stress-strain 

response did not. stabilize. This conclusion may be cycle and/or material 
model dependent. 

4. The results in Figure 36 show the cyclic hardening that occurs for-the 

one dimensional analysis during the first 15 loading cycles. 

5. — The inelastic strain components for cycles 1-15 are shown In Figure 37. 

As indicated, the creep strain per cycle remains essentially constant 
while the cyclic plasticity is decreasing. This Is consistent with the 
results from the 3D louver analysis and suggests that additional cycling 
would have produced a smaller amount of cyclic Inelestic strain tharv 
predicted for the 6th loading cycle. 

6. A comparison between the predicted response for cycle 15 and the results 
of the specinrten test Is shown in Figure 38. The prediction has a larger 
tensile mean stress and a smaller amount of cyclic inelestic stral i 
(defined as- the width of the loop across the strain axis). 

It appears that further modification to the material model (combined hardening 
and creep simulation) Is required to accurately reproduce the Hastelloy X 
thermomechani ca-1 response. 

4.2.6 Linear El astic Analysis 

Within the design process, rigorous non-linear analysis as described In Sec- 
tion 4.2.4 may not be alwa^ys practical. This Is particularly true where sever- 
al iterations are required to establish the final design geometry. A linear 
elastic analysis of the 3D louver model was conducted to establish the rela- 
tionship between the elastic response at the louver lip and the results pre- 
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(tlcttkl l).v Uh' nonlliHMr analysis. Tlio nosults wpt'o usinl to tdiMitify o-simpli- 
f Uxi (short uit') approach to (pial Huti vely detormino the actuaJ rosponso t row 
tho elastic response. Tiw' anaTysis used the same finite elauent- model and 
incremental t«npi>rature loading as the nonlinear analysis. The yield points in 
the stress strain curve description were set" at an arbitrarily hiqh value to 
prevent yield! nq durihq the cycle. — liUi— cxeep simulation was- not included in 
the elastic analysis. 

Predicted stress-strain results at the end of the lip are shown in Fiqure 3d. 
Aqain usinq the saiite letter desiqnation as the non-.linear analysis the 
response can be foUowed; 

A - Start of cycle 

C* - Minimiitt mechanical strain poi?it 

D - End of heat-up portion of cycle 
F - Maximiift mechani cal strain point 

A' - End of cycle 

A comparison of these results with tlw; predicted non-li-near analysis results 
indicates: 



Elasti c 

Non-Linear (6th C 

Max. Strain % 

+.0135 

-.100 

Min. Strain % 

-.3f?l 

-.448 

Strain Ranqe 

.3751 

. 348 

Mean Sti'ain % 

-.174 

-.374 

Mean Stress (MPA) 

- P45 

+ 138 

(KSI) 

(- 35.5) 

(+ 30.0) 


The cyclic strain ranqe p»'edicted with the elastic analysis is within 8% of 
the value from the fdh cycle of the non-linear analysis. While the peak stress 
and strain values are different between tlw? analyses, the shapes of the 
strai n-tem|.x'rature histories are similar. These results, shown in Fiqure 40, 
are due to tlie fact that the louver is primarily a theniially loaded structire 
and the variation in theniial <:train (*' .iT) ,..wMc!l produces the stresses, is 
indt'peruient of material stiffness. 

TIk? clastic strai n-tt?m|.H'rature history was used in the development of a quali- 
tative approximation of the non-linear results from the elastic calculations. 
Refw'inq to F iquro 41, the predicted mechanical strai n-ttxiiperature history and 
the monotonic stress-strain curves (Fiqure 15) are used to construct the 
stress-strain response diminq tiu' loadinq (heatinq) portion of the- first load 
cycle. This corresponds to the minimum strain point (0) in the Fiqure. Flastic 
un loadinq is assuiied to occur froin this point at an averaqe elastic modulus. 
Estimation of the reverse yield point is based on two- ass imp t ions. 

1. All of the mechanical strain qenerated durinq the heatinq portion of 
cycle is plastic strain. 



Thi' sizo of the yield surface durimi Bn1oad1nft.. f . m: .,... ea ch ti'inperaUiiro is 
doUwidned using an isotropic rule. 

Ttw rt^verse yield point is labeled 1; in figure 41. Continued unloading (cool- 
ing) takes place along a<» average i^vlastic modulus^ determined frimthe temper- 
atures in this portion of the cycle, to point f (the maxiimm strain). Reload- 
ing for the second cycle again usGs an average, el as tic modulus. The new yield 
point is detennined with the same a^*sumptions, however, the- equivalent plastic 
strain has been increased based on the revers^Kf yielding during thtr imloading. 
This- produces an atiditional amount of strain hardening which raises the new- 
yield point for the loading portion of the cycle, Cewparison of the cyclic 
response using this approach with the ?nd cycle response from the non-linear 
analysis is- also shown in the- fiiuire. Further cycling would produce a more 
elastic response (small inelastic strain range) due to the increased yield 
stress associated with the isotropic hardening assumption.. Strain rachetting 
cannot occur because the strain history is fixed by the elastic solution. Time 
dependent (creep) effects- have not been included.. This approach has produced 
good qualitative results, for the one dimensional stress state under strain 
controlled loading, relative to the more rigorous non-linear analysis. It may 
prove useful for future analysis of this typo where an approximation ot the 
response is acceptable. 
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Figure 12 MARC Finite Element Breakup of Combustor for Liner Specimen 
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Figure 14 Construction of Stress-Strain Input for Analysis 
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Figure 15 


Hastelloy X Stress-Strain Representation for MARC Analysis 




Figure 16A -Prediction of Cyclic Response at 76QOC (140QOF) 
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Figure 16B Prediction of Cyclic Response at 871QC (I 6 OOOF) 
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Figure 24 


Inelastic Strain Components versus Cycle Nunber 
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Figure 26 Louver 3-0imens1ona1 Model 
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Figure 29 6th Cycle Stness-Strai n 
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Figure 30 


Tubular Specimen for Thermomechanical Testing 
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Figure 33 Sensitivity of Specimen Response to Total Strain Range 


48 




AXIAL STRESS STRESS 


MPA 

(KSI) 



Figure 35 Comparison of Predicted 6th Cycle and Experimental Response 
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Figure 36 Cyclic Hardening Predicted by One Dimensional Analysis 
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Inelastic Strain Components From One Dimensional Analysis 
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5.0 TASK III CYCLIC LIFE PREDICTION AND SENSITIVITY ANALYSIS 

5.1 Introduction 

The cyclic thernloraechani cal loading experienced by the 
representative of the loading on many of the structural 
section of the gas turbi^ie engine and is one of the most 
for accurate durability prediction. The well 
tical results obtained for the louver component have been used 
the evaluation of several existing life 

mode for the louver has been considered as two mechanisms, crack 
the edge of the louver lip and axial crack propagati^in f 9 © f tbe 

11d toward the weld. Three crack initiation procedures, described in the cur- 
rent literature, were selected and applied to this pnoblm to determine their 
capability to predict the louver Tip crack initiation life. A linear ®^stic 
fracture mechanics analysis was- conducted and the results compared to-^the 
observed life growth data. Finally, a sensitivity analysis was conducted on 
each of the three crack initiation models. 

All of the louver lip information presented in the previous sectiot^ has been 
predicted at the integration point closest to edge of the lip. 
sis using an axisymmetric louver model to determine optimim element me^ 
density indicated that the strain gradient in the lip could be 
dieted by the 8 reduced integration point elements selected for the analysis. 
For the crack initiation life predictions, the stress and strain 
corner of the louverJIp were considered as the major driver for the initia- 
tion process. 

The extrapolation procedure used to determine these values 
last 3 lip elements and is shown in Figure 42. The 

layer of integration points was used with a extrapolation to deter- 

mine the surface strains. The two values of surface strain were then linear y 
extrapolated to the corner of the lip. Using this procedure, the 
between the corner node and the integration point in of total mechanical 

strain range is 7 %. A similar procedure for the inelastic strain 

(defined as^the tension going plastic strain from the ?n?ear^ 

a 17 5?!^ increase in the value at the corner relative to the nearest 

ti on point. In the following life prediction calculations the extrapolated 

corner values are used. 

5.2 Crack Initiation Analyses 

5.2.1 Strain Range Partitioning Method 

Strain Range Partitioning (Reference 5) is a well 
diction of^high temperature (creep/fatigue) crack 
is based on the assumption that a structural response contains time 
ent plastic and time dependent creep conponents. Under c^ic reversed load- 
ing, these component’s combine to produce four (4) possible types of damage 

modes: 
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pp " tensile plasticity reversed by coinpressi-ve plasticity 
pc - tensile plasticity reversed by compressive creep 
cp - tensile creep reversed, by compressive plasticity,, 
cc - tens41e creep reversed by compressive creep. 


A complex structural response is considered as containing as many as three (3) 
of these damage modes. The usefulness of this method has been, demonstrated bv 
many researchers .^on isothermal test specimens subjected to complex loading 
histories (Ref. f;. Life prediction by Strain Range Rartitionlng reoulres 

( A-cvs.. Life^ for the 4 generic damage- modes^ of the 
considered, Ln addition, the inelastic strain components gen- 
erated by the actual loading cycle in the component must be determined. For 
^e leuver. the predicted results, at the end. of the lip were used to define 
ttie magni tides of the damage modes. The thermomechanical response was assumed 
^ fP (^^sile plasticity revised by compressive plasticity) and 

*=’®vised by compressive creep) damage. This is cottsist- 
u}^? assumption that all of the time dependent response is occurinq at 
the higher -temperature when^the lip is in compression. 

The generic pp and pc life lines used in the life prediction were generated 
from a separate Pratt & Whitney Aircraft funded fatigue test program. The pro- 
gram consisted of uniaxial, strain controlled tests of fine grained Hastellov 
X material between 5380C (lOOOOP) and 9820C (I8OOOF) Fully reversed 
compressive hold time (pc) tests were included. The strain rate used 
for the pp tests (0.1 min-1) was considered to be fast enough. to prevent any 
creep or relaxation during these tests. The majority of the pp tests were con- 
ducted at 8710C (1600“F), however, the limited number of specimens at the 
Other temperatures did not indicate a strong temperature dependence on fatigue 

lilalp^5;,f. PP w®s considered as a 

single data set and a least squares fit used to define the life line. In the 

compressive hold time (pc) testing, the same transient strain rate (0.1 
min 1) was used to reverse the strain and a 1 minute strain controlled hold 
time imposed at the maximum compressive strain point. After reviewing th2 
data. It w^ determined that the inelastic strain generated in these tests 

thS^ tl^ f damage ( AW Axm =0.5 - 0.75). This required 

interaction damage rule and the pp life relationship be used to 

7finor * /^ number of pc tests were conducted at 

760 C (1400 F) and 871 OC (1600°P)* The results showed a stronger tem- 
perature dependent than in the pp tests, however, all data points were 
hT to establish the pc life relationship. Tte prand jc 

life lines determined from these tests are shown in Figure 43. Their equations 

dr Qm ” 


Npp = 




-1/.5I 


5.1 


The Initial life prediction with the SRP method used, the results from the 6th 
loading cycle of the non-linear analysis. The pp and pc components were deter- 
mined from the amount of compression going plastic and creep strain as pre- 
dicted by MftRC. The inelastic strain components at the element integration 
point are Aepp » 0.025^ and A.e pc = 0.06*. Extrapolating, the strains to the 
corner of the louver lip and using the above life relationships with the 
interaction damage rule produces: 


1 

Fpp 

I'pc 



Npp 

Npo 


Acpp = 

.o:94'; 

Aepe 

= .0705 

I'pp = 


I'pc 

= .700 

Npp = 

lOOOO 

Npe 

= 7850 


8500 oyd 

os vs. 1000 

1 :50 cycio: 

ohsorvcil life 



This overprediction of the crack initiation life may be associated with un- 
certainty in the definition of the fatigue life relationships (particularly 
the pc component). Additional fatigue data could potentially reduce the pre- 
dicted life, however, it appears that the SRP method would still overpredict 
the cracking by at least a factor of 2. 

5.2.2 Pratt & Whitney Aircraft - Commercial Products Division Combustor Life 
Prediction Method: 

Development of this method is fully described in Reference 7. It is concept- 
ually the sane as the Strain Range Partitioning method in that the material 
cyclic damage- is assumed to be parti oned into separate "types". For a combus- 
tor loading cycle the PWA system as-sunes that the damage contains only 
pp- (compressi ve plasticity reversed by tensile plasticity) and pc-(comp‘.'*essi ve 
creep reversed by tensile plasticity) components. The inelastic strain <"ange - 
life relationships used are similar in form to the "ductility normalized 
relationships proposed for SRP. They are: 


Npp 


/ Arpp 
.75 \ Dp 


■1 (),(. 


Npi- 


yi.75 \ Dp 


-l/O.S 
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where: Npp 

Atpp 

Dp 

Npc 

^epc 

A linear damage 


- fatigue life produced by the pp component 

- pp component of inelastic strain range. 

- temperature dependent ttDsile ductility 

- fatigue life produced by the pc component 

- pc component of inelastic strain range 

rule is used to predict total 1 ife. 


1 1 ^ I 

Ntotal Nw Npc 


Prediction of the inelastic strain range components (^^pp, ^^pc) for this 
method is based on the total mechanical cyclic strain range determined from a 
linear elastic analysis and two algorithms: 


Acpp = At\.| 




where: ^ ^ TR 


total strain range predicted from linear elastic analysis 


Acgi 


endurance limit strain range associated with a fatigue 
life of 10^ cycles 


E - Young's modulus 

Ki - temperature dependent constant used in stress- strain 

curve representation ( e = K^a 

- time and temperature dependent constant developed from 
creep data for prediction of relaxation response 
{ f-c - Bfo "*) 

Prediction of the louver lip crack initiation life was based on the following 
val ues4 

= 0.401% - extrapolated total elastic hoop strain range at 
corner of lip 

On * 38% - tensile ductility associated with maximum cycle 

^ temperature (9820C (1800OP)) 
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Kl= 6,8 X 10"7, 3,1 

Bt =5.7 X- 10“5 - Bas4d-on strain hold time at 9820C (iaOO°F) 
E = 18.4 x -103 KSI 


^^el * .18% - endurance 1 imft strain 
temperature 9820C (1800QF) 


range at maximum cycle 


The resulting quantities are: 

= .221%; ^fpc » .068% 
Npp = 3300 ; Npc = 3600 


Ntotal = 1720 cycles 

This method predicts a life closer to the observed component life, however, 
the predicted inelastic strain range { Aepp + Aepr) is nearly 2.9^ times 
as large as the inelastic strain observed in both the non-linear analysis and 
thermomechanical uniaxial speciment test. Reducing the inelastic strain range 
to the observed value increases the calculated life to approximately 8000 
cycles, cons istnt with the SRP prediction. In actual design practice calcula- 
tion using this method would be used with experimental engine and field ser- 
vice data to determine an overall service life of the component. 

Both of the above methods, based on the prediction of inelastic strain, over- 
predicted the louver cracking life. The fact that each is calibrated using 
isothermal fatigue and creep/fatigue data suggests that thermoroechanical 

(varying temperature) cycling may produce fatigue damage at a aster rate than 
an isothermal cycle. 


5.2.3 Continuous Damage Method: 

The third crack initiation life prediction for the louver lip used the Conti n- 
uous Damage Approach proposed by J.L. Chaboche (References 8 and 9). This 
method was selected for evaluation because it is sufficiently different from 
the previous two methods (P&WA-CPD Combustor Life Prediction Method and Strain 
Range Partitioning). The approach is similar in concept to the Ductility 
Exhaustion method used in the prediction of turbine blade creep/f atique life 
at Pratt & Whitney Aircraft (Reference 10). 

The Continuous Damage Method is a stress based approach which assumes a damage 
internal variable that defines the strength of a material.. The prediction of 
life requires solution of differential damage equations which determine the 
evolution of this damage. The damage (D) is assumed to be 0 in the initial 
und^aged state and 1 at failure. For high temperature loading the damage con- 
sists of creep (time dependent) and fatigue (time independent) components. 




The creep damage component follows the fom of K'acharvov (Ref. 11) and Rabot 45 ov 
(Ref. 12) and is written: 


ill) 


/ lf>l \' -l; 
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where: D = damage variable 

A,rji_ = temperature -dependent coefficients 

t- a time 

The fatigue damage, process is defined as: 
dD = 1 1 (1 1 
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where: 


a,.b,Mo,^ 

a 

temperature dependent coefficients 

0 m 

a 

maximum stress in the cycle 

0 

= 

mean stress in the cycle 

ou 

= 

ultimate tensile strength 

01 

s 

fatigue limit under completely reversed loading 

0 1* 

S 

fatigue limit for non-zei^o mean stress 

N 

a 

cycles 


This expression describes, for rapid cycle fatigue (no creep), the relation- 
ship between maximum stress and life including mean stress effects. Combining 
the two damage expressions results in the interaction equation: 

( <1 l >> ^ 1 1 (1 1 ))^^' 1 ^^ ) 

|lk+l)Nc (^?+l)N,, (l cvi(l ) 


ill) = 
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where 

lie 

Nf 


= total life 

1 

= life in pure creep = “* = 

N (;’' 


<e f ati gue 



(k+ 1) 


rM 


= ( 0 + 1)0 


dt 


om a 


Mo[l - b(-)l 
- all 




h ' 
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■or the combustor determined from the high 

the creep portion of ^ constants asso- 

stress-short time creep were obtained from the 

ciated with the fatigue portion of the Aepp vs. Npp 

t-rteJlirtfe c=f. « were m^de for the 

louver life prediction with this method, they are* 


1 . A ..ean .tress has 

S 2 trst ' re : s ’‘’ v ^“: es '^<^^38^ i .' A * (^ f ) 

simolifies the fatigie equation allowing. 


h = f 

ttV* = ( 

M(ti)= 1 


no 
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Z, -The creep (relaxation) damage 1& assumed to occur only between 871°C 
(1600°F) and 954®C (1750°F) In- the louver cycle. Based on the 
creep testir^, t"he constants A, r,. K are temperature independent.. No 

difference is assumed to exist between -tens i on vs. compression and 

creep vs., relaxation damage.. - 

3. For the case of a variab-le temperature history (thermomechanical 
fatigue), Chaboche proposes the idea of an effective temperature which 
is- used to evaluate the constants in the fatigue equation. The effec- 
tive temperature is determined by considerijig the themoraechanical 
fatigue cycle on a -non^d-imeTiSionalized stress basis which— results in 
the equation: 


lid'll 



i3(T)<dl2> 

Mo T, IN, HI h M,r,d') 



1 1 ~ rediu'ed moan stros.s 

mi 


r = toinpcruturc 

I'* = olToctiw tomporatmv 

'^1-(|M MM,. D = isolluMinal latifjiio litV 


Nj. (|’||. ||, T'*=) = latiiiuc lit'o iiiulor varying 
tomporatiiiv 


For a given stress -temperature history and isothermal fatigue life relation- 
ships the integral can be evaluated. The effective temperature (T*) is then 
determined from the maximim stress amplitude (Pn) and the value of Np 
(Pi I* Il» T*). This approach has been demonstrated by Chaboche (Reference 
9) on a high strength, fairly rate insensitive material (IN 100). Following 
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tMs. proceduce for the HasteUoy X Isothermal fatigue tests produeed. a- reduced 
tatigue cucA/e (I? vS N) which was temperature Independent wh-lch- makes the 
concept of an effective tempeciat ure 1 nappropri a-te for th& Hast&loy X thermo- 
mechanical cycle. Furthermore, the stress levels obtained the thermomeGhan- 
ical cycle are smaller than those obser-ved in the rapid cycle (pp) fatigue 
tests. This is due to the high rate sensitiviity of the- mater i at at these tem- 
peratures, Those results results, showft in Figure 44, suggest that the majori- 
ty of- the total dJinage In the louver cycle is associ-ated- with creep and not. 
time independent fatigue. These observations raise sane doubt as to the appro- 
priation of this method for the thermomechanl cal fatigue of a material nke 
Hasten oy X. 


In the actual life calculation for the louver Vip,. several fatigue damage com- 
ponents were* considered to determine the sensitivity, of the caVcal ation,. The 
predicted 6th. cycle stress-Strain response at the end of the louver lip was 
used to establish the stress levels for the damage epoations (see Figure 45). 
The creep damage is assumed to occur between-- 151, 7 MPA (-22 KSI) and -z/.b 
MPA (-4 KSI) within 1 minute, an exponential die out of stress with time is 
used. The stress amplitude is chosen, as 358.5 MPA (52 KSI) (the maximun ten- 
sile stress) and 510.2 MPA (74 KSI) (the 

fatigue damage is assumed to occur at either 538PC (lOOOdp) or 760 C 
(1-400°F). A summary of the material constants and life predictions are shown 

bel ow: 


Hasten oy X Constants for Continuous Damage Approach 


Temperature (°F) 

(KSI) 

1000 

86 

1400 

57 

A = 7.8 X 107... 

r = 0.23 

Stress Afflolitude 

F atigue Life 
(cycles ) 

358.5MPA @ 5380C 
(52 ksi 0 lOOOOF) 

128300 

510.2MPA 0 5380C 
(74 ksi 0 lOOOOF) 

600 

358.5M5A @ 760OC 
(52 ksi 0 1400OF) 

745 


°1 (KSI) 

(3 

Mo 

42 

8.92 

219 

32 

6.51 

242 

K = 3.51 

a = 1 

b = 0 

Creep Life 
(cycles) 

Total Life 
(cycles) 

Actual Life 
(cycles) 

515 

495 

1000^1250 

515 

250 

1000-1250 

515 

250 

1000-1250 


In all calculations the creep component was the limiting life component. Addi- 
tional variations In the stress levels during relaxation did not signyicantly 
Increase the predicted creep life. The conservative prediction of the louver 
lip cracking life suggests that the assumptions of the stress amplitude pro- 
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ducing fatigue damage and damage- associated with compressive, relaxation- may 
not be appropriate for Hastelloy )i at elevated temperature. This Is further 
demonstrated when, the Contdnuous Damage Method was used to predict the life of 
the Isothermal strain contro1l.ed hold tiros- tests. These are the. same creep 
fatigue tests used to generate the Aepc vs, N. curve for the Strain Range 
PartUlonIrvg Method. For all tests (5)» the Continuous Damage Method, predicted 
life was approximately IS% of the actual life. This concluislon Is consistent 
with the general SRP observation that a cp (tensile creep- compressi ve plasti- 
city) is more damaging than a pc (tensile plasticity- compressive creep) type 
of cycle. 

5, 3 Crack Propagatl on Anal ysi s;. 

A ser-les of analyses were conducted- using the linear elastic fracture mechan^ — 
ics method to predict the subsequent growth, of fatigue cracks- from the edge of 
the louver lip,. Existing Hastelloy X sheet crack grow.th. data was. used to 
establ ish. the crack growth. Taw for the analyses. This data was obtained from 
isothermal, center cracked. panel, load controlled (R » 0) tests which had been 
conducted at temperatures between 427<^(T" (800°F) and 760°C (1400°F)., 

Only the 760O(r (140QOF) data waS used in the analysis, since It represent- 
ed the fastest available growth rate and would provide the most- conservative 
prediction. This assumption Is consistent with the louver lip stress- strain 
response since the highest tensile stresses occur during the therrncmechanlcal 
cycle at temperatures between 760°C (1400°F) and 504^0 (940®F). The 

louver lip was modeled as a semi-infinite plate with a through thickness edge 
crack. The prediction was based on an Influence or weight function method 
described in References 13 and 14. This method allows for very general stress 
(strain) gradients to be simulated. The crack driving force (AK) Is computed 
from the stress (strain) distribution in the uncracked part, along the plane 
of crack propagation. The effects of the cyclic stress range on the predicted 
growth rate was modeled by means of the R-ratlo. The general equation for the 
crack growth rate used In the analysis_lsi — 

[_ m 

fL] 

1 Rj 


where: c,m » temperature dependent constants 

R = R-ratlo ( °m1n/ ‘^max) 

A-K = stress intensity factor 

The stress range ( ao) and R-ratlo In the analyses were determined from the 
6th cycle of the louver non-linear analysis. The distribution of these quanti- 
ties along the louver lip is shown In Figure 46. The results of three (3) cal- 
culations are compared to the measured rig crack growth data In Figure 47. In 
all calculations, an Initial crack length of .0254 cm (.010 In) was assimed to 
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exist at 1050 cycles. Curves A and B ar-c the predictiORS as^^ing R = 
sidering the. e^l tire stress rat^ge as. te-nsUe loading^. The difference |j- 
two curves is associated- with two representations, of the crack growth law 
shown-in the Figure. Curve A used the UUnear representation while B-resulted 
from the extrapolation, of the low A K values. This modi-f icati on to the crack 
growth law was investigated because the steeper portion, of the law (large ak) 
was associated with net section yielding ahead of the the center, 

cracked tests. Curve C s the result of including the predicted-R ratio dis- 
tribution along the lip in the crack growth, law. The negative R-raiio reduces 
the growth rate (da/dn) for a given.- AJ< resulting in much slower, response 
than^the R » 0 calcul ati ons . All 3 calcalations predUt an initial crack 
growth rate which is slower than the rig data. Curve A canes the closest to 
predicting the cyclic life, but the growth rate at crack l-engths greater th^ 
.0508 cm. (.060 in) is much steeper than the data. The effects of time depend- 
e.nt plasticity and the varying temperature cycle have not been- included_ in 
these predictions.. Under these conditions^ the local stress_ at the crack tip, 
affected by closure and stress relaxation, may accelerate the predicted crack 
growth. rate. 

5.4 Sensiti vity Analysis 

A sensitivity analysis was conducted on the three crack initiation lite pre- 
diction methods considered in Section 5.2. Variation in total mechanical 
strain range was selected as the key parameter for assessing the sensitivity 
of each of the life prediction methods. An increase in the mechanical strain 
range in the louver geometry is- possible by perturbations of the temperature 
distribution and magnitude along the louver lip. The axi symmetric louver model 
oreviously used to establish the finite element mesh density along the lip, 
was used to predict the increase in circumferential (hoop) strain with these 
changes. Two variations in temperature distribution and a lOQOF increase in 
maximum temperature were considered. Linear elastic analyses 
using the 10 second transient temperature distribution as the baseline. This 
coincides with the point of maximum temperature gradient (strain) m tpe 
louver cycle. It was assumed that changes in the predicted strain at this 
point are directly related to changes in the mechanical strain range, 'be 
results of the analyses, (Figure 48) indicate that, as more of the tOTf^rature 
gradient is concentrated toward the edge of the lip, (profiles A and B) , the 
predicted hoop strain increases relative to the baseline value. An increase in 
temperature of lOQOF at the lip edge further increases the strain relative 
to the baseline. Although it is the more severe condition, profile B is least 
likely to occur due to heat conduction effects. Assum.ing a possible error of 
lOQOF in the temperature measurement during the test and the "A" temperature 
profile, an increase in the total mechanical strain range at the edge of the 
louver lip of 20X seems reasonable for the sensitivity analyses. 

The sensitivity of the three life prediction methods to a change in total 
mechanical strain range is shown in Figure 49. These calculations reflect only 
a change in the strain range, any change in temperature dependent parameters 
is neglected. Each method will be discussed separately. 
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PWMIBD Combustor L-ite Pre diction Method; 

This- method has a moderate DredicU^^ 

a t,f s?r\Tr*4n;^.7e&^^ the two 

JlJSritiSis iscuased fn Section 6..2.2, Increaaes >>2 Jf 

cliSge U associated with Aapp damage. A sunmary of the two calcu.1 attons 

shown belcw. 


Basel ine 


120?^ of Baseline 


total 

inelastic 


pp 

A e pc 
Npp 
Npc 
N total 


0.405S 

0.289% 

0 . 221 % 

0 - 068 % 

3299 

3592 

1710 


total 

Inelastic 

f, pp 

pc 

Npp. 

Npc 

Ntotal 


0.481% 

0. 373% 

0.301% 

0.072% 

1964 

3336 

1240 


Strain Range Partitioning ; 

Janae^ The results of the uniaxial specimen testing shown in demon- 

s“?t« this effort. These ^““1“ mphaslz^^ 

M?slra1nTanrJo™A^tf^p"^ * calculations 

are shown below. 


Basel ine 


120% of Baseline 


total 
inel astic 


Nop 

Npc 

Ntotal 


0.372% 

0 . 100 % 

0.029% 

0.071% 

10636 

7847 

8500 


to 

in 

a! pp 

pc 

Npp 

Npc 

Ntotal 


total 

inelastic 

PP 

pc 


0.446% 

0.174% 

0.051% 

0.123% 

4700 

2640 

3000 


Continuous Damage Method : 

*dfcUd’ikl,’:aMnrthe“pr%dl"rt%'^^^^^^^^ fP'«1ve to the data. 
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Figure 42 


Extrapolation of Louver Lip Strain 
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Figure 43 General Hastelloy X PP-PC Life Relationships for Strain 
Range Partitioning 
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Figure 45 


Stable Stress Response for Continuous Damage Method 




Figure 47 Canparison of Test Data and Crack Propagation Analyses 







F igut'e 49 


Sensitivity of Three Crack Initiation Life Prediction Methods 



6.a DISCUSSION OF RESULTS AND-CONCLUSIONS 

This- program has addressed most of- the disciplines associated with, the design 
and analysis of gas turbine engine hot section structural components. The 
specific analyses have been performed- on a well controlled simulated combustor 
liner test specimen.- These results have identified a number of areas where 
further research and development would contribute toward the improved durabil- 
ity of these components. General conclusions and observations concerning this 
work are: 

1_ Wot section structtiral components are S 4 jbjected to complex thermomechan- 
ical loadings. A well controlled source of benchmark data is a useful 
tool for the development and calibration of analytical and life -predic- 

tion PK)dels. 

2. The three dimensional finite difference heat transfer solution conducted 
for this program accurately predicted the louver steady state and tran- 
sient, temperature distribution. 

3 ...... The temperature predictions are highly dependent on the accuracy of the 

thermal boundary conditions. The greatest uncertainty in this analysis 
was in the calculation of the louver heat transfer coefficients.. 

4. The three dimensional spatial and temporal interpolation between the heat 
transfer and finite element structural meshes required a large amount of 
effort. The interpolation scheme used had difficulty in consistently 

selecting thermal nodes to determine a temperature for a structural inte- 
gration point. A number of the temperatures were manually adjusted to 

agree with the temperature solution. Part of this problem may have been 

associated with the louver dimensions and the small included angle in the 

model. Potential improvements in this area would reduce the overall man- 
power to complete a thermal-structural analysis. One area suggested for 
investigation is the use of a finite element heat transfer solution. Use 
of a consistent element type and mesh would then eliminate the need for 
interpolation. 

5 The combined (isotropic-kinematic) hardening rule can be used to model 
the isothermal Hastelloy X strain controlled cyclic response at a given 
strain rate^ 

6. A constant creep strain rate model developed from short time-high stress 
level creep tests accurately predicted Hastelloy X isothermal stress 
relaxation data. 

7. Results of the nonlinear louver structural analysis indicated that the 
time independent plasticity model and the creep simulation did not accur- 
ately predict the cyclic thermomechanical response at the louver failure 
location. Tests of a uniaxial strain controlled specimen run with the 
same mechanical strain- temperature history as predicted at the failure 
location indicate that the stress-strain response stabilized within the 
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first few- cycles. Simulation of this experiment with the Hattelloy X 
model discussed above showed continued cyclic hardening (increasing peak 
tensile stress and reduced inelastic strain range) after many cycles. 
Potential modification to the material model, including a multi-yield- 
surface concept V non-1 inear hardening or use of one of the rate dependent 
(unified) theories currently under development, may be required to 
improve the prediction for the varying temperature loading condition. 
Determination of correct thermomechanical response is important in the 
life prediction, and in determining the. extent (number of cyoles) of noa 
linear analysis required for a str4jctural component. 

8. The three (3) h-igh temperature (creep/fatigue) life prediction methods 
considered are the- Strain Range Partitioning,. Pratt & Whitney A.ircraft - 
Commercial- Products Division Combustor Life Prediction, and Continuous 
Damage methods. All assume that time independent plastic and time depend- 
ent creep damage mechanisms are present at elevated temperature, l-so- 
thenndl fatigue and creep rupture teste are used to define the material 
life relationships. 

The Strain Range Partitioning and PWA-CPD methods assume the existence- of 
generic types of fully reversed damage cycles composed of canbinati-ons of 
the plastic and creep mechanisms. For purposes of this analysis, the com- 
bustor louver lip response contained only the- pp(tensile plasticity 
reversed by compressive plasticity) and pc(tensile plasticity reversed by 
compressive creep) damage cycles. 

The Strain Range Partitioning method overpredicted the louver cracking 
life (8500 cycles vs. 1000 cycles). Part of this discrepancy may be asso- 
ciated with uncertainty in the definition of the generic pp and pc 
fatigue life curves. Better definition of the curves may reduce the pre- 
dicted life (and improve the correlation), however, it appears that the 
SRP method will overpredict the cracking life by at least a factor of 2. 

The Pratt & Whitney Aircraft - Commercial Products Division method also 
overpredicted the louver cracking life (1700 cycles vs. 1000). This 
improved correlation (relative to SRP) is due, in part, to the fact that 
the inelastic strain range predicted by this method is larger than the 
observed inelastic strain in the louver.. Equating the inelastic strain 
value results in a predicted life of 8000 cycles, which is similar to the 
SRP calculation. In actual design practice, this method is used with 
experimental and field service data to assess the overall service life of 
the component. These results suggest that a thermomechanical fatigue 
cycle may produce damage at a faster rate than a comparable isothermal 
cycle. 

The Continuous Damage Method predicts life based on the growth of a dam- 
age parameter (D) from the initial undamaged state (D “ 0) to failure (D 
= 1). The method uses a non-linear accumulation of damage produced by the 
plastic (fatigue) and creep mechanisms. This method is a stress (rather 
than strain) based approach and has the general formulation to consider 
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;t;fse:Sti:/rri.a1l‘1nrt1af^cV,^ 

predicted, the louver Up cracking life by The 

riiiAtinns the creeo component was the controlling life component, ine 

assumption that compressive stress ^®^axation (observed in th^ 

cycle) is as damaging as tensile creep may be- partially responsible tor 

the conservative prediction. 

Q- rai<»iii Afi nn of fatioue' crack growth using the linear elastic fracture 
mechanics approach, may not be appropriate for Hastelloy X thermomechan^ 
ical loading. Analyses w6re conducted ® 

from isothermal (7600C (1400OF)), 

rracked oanel testing. Stress values were obtained from the 6th eye e 
nonlinear analysis. The use of load-controlled data and predicted stress 
Tanaef was ass^umed to be valid tor the small crack lengths considered. 

Two calculations were performed,. 1) assuming that the '^maxl'^dis- 

wAc ton<iile loading {R=0), 2) using the global R ratio ( min/ max) dis 

tribution determined from the louver analysis. Both 

tial growth rates considerably slower than the observed values. 

Inelastic and time dependent effects were neglected, in fbese calcuTa- 
tionf The Targe region of inelastic response, relative to the small in - 
tial crack sizes, suggests that % 

mav be a better correlating parameter. In addition, stress 

??Jn obse™d in Hastelloy X suggests that the local R ratio near the 

crack tip is more important than the global R ratio effect. 
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